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Abstract— This article proposes an algorithm for generation of N-dimensional rotation matrix R, N=m+n, 


m=2?, n=24, p,q € [2, 4, 8] which rotates given N-dimensional vector X in the direction of coordinate axis 


x; Algorithm uses block diagonal matrix, composed by Transpositions matrices. As practical realization 


article gives Matlab code of functions, which creates Householder and Transpositions matrices and V 


matrix for given n-dimensional vector X. 


Keywords— Applied Mathematics, Mathematics of computing, Mathematical analysis, Numerical 


analysis, Computations on matrices. 


I. INTRODUCTION 


This paper proposes an algorithm for creating an n- 
dimensional rotation matrix V, which rotates a given 
vector X in the direction of the coordinate axis x; using 
block diagonal Transpositions (Trs) matrices. Matrices that 
rotate a given vector X in the direction of a given 
coordinate axis are used in calculations such as QR 
decomposition and in obtaining a rotation matrix 
according to the NRMG algorithm [2]. 


Il. ROTATION OF N-DIMENSIONAL VECTOR X 
TO THE DIRECTION OF AXIS Xi 


The rotation of an N-dimensional vector X to the direction 
of one of the coordinate axes is a basic operation of the 
NRMG algorithm for obtaining an n-dimensional rotation 
matrix [2] and for QR decomposition [9] of a matrix A 
into a product A = QR of an orthogonal matrix Q and an 
upper triangular matrix R. 


Definition 1: We will call a V matrix or Vx matrix an 
orthogonal matrix that rotates a given vector X in the 
direction of the x; coordinate axis. 


To obtain a V matrix for dense vectors, the method using 
Householder reflections is considered the most effective 
[7],[13], while for sparse vectors the method using Givens 
rotations may be more effective. However, for vectors of 
dimensions 2, 4 and 8, the method using a Transpositions 
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matrix is the most effective [1]. The following is a 
description of these three methods for obtaining an 
orthogonal matrix V that rotates a given vector X to the 
direction of given coordinate axis (e.g. axis X1). 


2.1 OBTAINING AN V MATRIX USING 
HOUSEHOLDER REFLECTION 


According to the transformation of Householder [13] if X 
and Y are vectors with the same norm there exists an 
orthogonal symmetric matrix P, which rotates vector X to 
the direction of vector Y such that 


Y = PX where P=I-W.W! 
and W =(X - Y)/|X-—Y| 


(1) 
In particular, if Y is unit vector Y=[1,0, ... 0], we get 
Y =PX where P=I-W.W! 
and W = (X - Y)/|X -Y| 
X= xxl e 


then matrix P will rotates vector X to the direction of the 
coordinate axis x1 Matrix P is matrix of reflection (not a 
rotation) because det( P) = -1, (which gives the name to 
method) 


The following example shows the Matlab code of the 
function Householder(X,k) that creates and returns 
orthogonal matrix V which rotates given vector X to the 
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direction of the coordinate axis x, using Householder 
transformation. 


function V= Householder (X, k) 
m=length(X); 
k=fix(k); 
v=X; 
v(k)=X(k)- norm(X); 
s=norm(v); 

if s~=0 
w=v/s; 

V=eye(m)-2* w*w'; 
else 
V= eye(m); 


end 


Example 1 Code of Matlab function, which returns 
orthogonal matrix V, that rotates vector X to the direction 
of the coordinate axis xx using Householder 
transformation 


As it can be seen, function Householder(x,k) performs two 
normalization of vectors and calculates elements of matrix 
w.w!, which needs n? floating points multiplications. 


2.2 OBTAINING AN V MATRIX USING GIVENS 
ROTATIONS 


Rotation of given vector X in the direction of one of 
coordinate axes (e.g. axis xl) can be performed by 
subsequent multiplications by Givens matrices [2] as 
follows: 


1 
Xmag = [[G(k,k+1,0,).X=My.X= 
k=N-1 
[rx,0,...0,0]" 

(3) 
where G(k, k+1, 6k) is Givens matrix, which rotates vector 
in coordinate plane (k, k+1) by angle Ok and rx is the norm 
of vector X. 


Givens matrices G(k, k+1, 0k), k = N-1, n-2, ... 1 are 
defined as follows [4],[5]: 


1 0 >~ 0 0 «+ oO 0 
0 1 = 0 0 > 0O 0 
Gk k+1,6)=|9 0 Ck -Sk o o0 
0 S, C o o0 
o 0 o oO 1 0 
o o0 o oO o 1 

(4) 
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where 9k is the angle of rotation and coefficients 
Ck=cos(6k) and Sk=sin(@k) appears at the intersections of 
k-th and k+1 -th rows and columns. 


The target of multiplication by Givens matrix 
G(k, k+1, 6k) is to set to zero coordinate Xk 4+1 - Itis easy 


to find that equation Xk +1 =() and equations (5) are 


satisfied simultaneously when sin(@k) and  cos(6k) are 
calculated using formulas (5) [2]: 


n X 
sin(0) =- 7, 

VXkK tXķk=1 

= Xk ae eee 
cos(8,) = if Xx +Xk=] >0 


Xe + X 


sin(0,) =0, cos(®,) =1 if x2 +x; =0 


(5) 
Thus searched matrix Vx can be calculated as 
multiplication of Givens matrices 


G(N -1,N—2,0,_,), GN -2,N -3, 0N 5)... G(1,2, 0; ) 


as follows: 


N-1 
Vx = [IG(N-k,N-k+1 0, ,) 
k=l (6) 


where coefficients of rotation sin(@k) and  cos(@k) are 
calculated using formulas (5). Calculation of angles of 
two-dimensional rotations 0k really is not needed. If 


x? F x? +41 =0 then corresponding Givens matrix is 
equal to Identity matrix G(k, k+1, 0k) = I — no rotation. 


2.3 OBTAINING AN V MATRIX USING 
TRANSPOSITIONS_MATRIX 


Definition 2 Transpositions matrix (Tr matrix) [1] is 
square matrix with dimension n=2™, m e N, which rows 
are permutations of elements xx, k €e [0, n-1] of given n- 
dimensional vector X and the values of elements Trp, q of 
matrix are obtained from the elements of given vector as 
follows: 


Troa F X peg p.q € [0,1,..., n =| (7) 


where © denotes operation ‘bitwise exclusive or’ (XOR) 
[15],[16],[17]. 


Definition 3 Transpositions matrix with mutually 
orthogonal rows (Trs matrix) [1],[2] is square matrix with 
dimension n=2,4 or 8, which is obtained as Hadamard 
product [10], (denoted by ©) of Tr matrix and n- 
dimensional Hadamard matrix whose rows (except the first 
one) are rearranged relative to the rows of Sylvester- 
Hadamard matrix [10] in order R=[1, r2, ... mJ’, 1, ..., 
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TnE[2,n] for which the rows of the resulting Trs matrix are 
mutually orthogonal. 


By definition, Trs(X) is a matrix whose rows are mutually 
orthogonal and contain transpositions of a given vector X, 
while the first row contains the values of X without 
inversions. Therefore, if ||X||=1 then the multiplication of 
Trs(X/|X||) by X gives unit vector to the direction of the 
coordinate axis x;. The preparation of an orthogonal Trs 
matrix from a Tr matrix and a Sylvester-Hadamard matrix 
for a given vector X is given by the formula: 


Trs(X)=Ts(X)o H(P) 
X=X/|x|| 8) 
where: 


e o denotes operation Hadamard product of 
matrices [13], 

e H(P) is n-dimensional Hadamard matrix, which 
rows are interchanged against the Sylvester- 
Hadamard matrix in given order P=[1, r2, ... taJ", 
T2, ..., f€[2,n] for which the rows of the resulting 
Trs matrix are mutually orthogonal. 

e X is the vector from which the elements of Tr 


matrix are derived. 


The formula below shows obtaining of 4-dimensional Trs 
matrix for a vector X=[x1,x2, x3, x4]? and for H(P) rows 
ordering P=[1 2 4 3]". 


Trs(X) = H(P)oTr(X)= 


ENTIT] Xi X2 X3 X4 
1-1 1 =l |X. x; X4 x3 
kee] j X3 X4 XX, X 
1 1 -1 -1) (x4 x3 X% x 


Xı X2 X3 X4 
Xə? —Xy X4 =X3 (9) 
X3 7X4 =X] X2 
X4 X3 =X = Xi 


It is important to note, that the ordering P of Hadamard 
matrix’s rows (against the Sylvester-Hadamard matrix) 
does not depend on the vector X but only on its dimension. 


An algorithm for creating Trs matrices for an arbitrary 
vector X for dimensions 2, 4 and 8 is proposed in [1]. 


The following example shows the code of a Matlab 
function that creates and returns Trs matrix which rotates 
given vector X of dimension n = 2, 4, or 8 in the direction 
of coordinate axis xı. 


H8=[l1111111;1-11-11-11-1;1-1-11-11 1-1; 
11-1-111-1-1;1-11-1-11-1L1;11-1-1-1-111;1 
-l-L11-1-11;1111-1-1-1-1); 


N=length(X); V=zeros(N); 


if N==1 V=X; 
else 
if N==8 H=H8; 


else if N==4 H=H4; 
else if N==2 H=H2; 
else error(‘length of X must be 2,4 or 8'); 
end; end; end; 
for col=1:N 
step=col-1; %Step of displacements 
for row=1:N % Loop for rows 
new_row =bitxor(row-1, step)+1; 
if H(new_row,col)==1 %set Trs value 
V(new_row,col)=X(row); 
else 
V(new_row,col)=-X(row); 


end; end; end; end 


function V = TrsMatrix(X) 
H2=[1 1; 1 -1];H4=[1 1 1 1; 1 -1 1 -1; 1 -1 -1 1; 1 1-1-1]; 
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Example 2 Matlab function, that creates and returns a Trs 
matrix for a given vector X of size n = 1, 2, 4, or 8. 


As it can be seen, function TrsMatrix(X) obtains elements 
of matrix V, which needs n? bitwise multiplications of the 
indices of the vector’s elements. 


II. COMPARISON OF THE COMPUTATIONAL 
EFFICIENCY OF THE METHODS FOR 
CREATING AN V MATRIX 


We will compare the computational efficiency of the 
methods for creating a matrix that rotates an n-dimensional 
vector in the direction of the x1 coordinate axis (V matrix) 
using Householder 
Transpositions _matrix. 


3.1 THE COMPUTATIONAL EFFICIENCY OF THE 
METHOD THAT USES HOUSEHOLDER 
REFLECTION 


Reflection and using 


Rotation of given vector X in the direction of one of 
coordinate axes (e.g. axis Xı) can be performed by 
subsequent multiplications by Givens matrices [2] as 
follows: 


As can be seen from formula (2), to obtain a matrix P 
(which for us is the required matrix V) the vector Y must 
be a unit vector Y = [1,0,..., OJ’. The calculation of the 
matrix R includes the calculation of the norm of the vector 
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X, the norm of the vector W = X / || X || - Y and the 
calculation of the elements of the matrix W.WT, which 
needs n? multiplications. Other computational operations 
due to their lower weight can be ignored. So, for the 
volume Cy of computational operations when using 
Householder Reflection we can write the following 
estimate: 


Cu = 2 norms + n? multiplications ............ (10) 


3.2 THE COMPUTATIONAL EFFICIENCY OF THE 
METHOD THAT USES TRANSPOSITIONS 
MATRIX 


As can be seen from formula (6), obtaining a 
Transpositions matrix Trs (X) for a given vector X with 
dimension 2, 4 or 8 includes the following operations: 


e normalization of the vector X 


e obtaining a Transpositions matrix Tr (x) for vector 
X =X/|X||. 


e element-by-element multiplication (Hadamard 
product) of the matrix Tr (x) by a matrix H(P) 
having the same dimension. 


As can be seen from the definition of the Trs matrix (2,3), 
its production includes n? bitwise multiplications of the 
indices of the vector elements, and (since the matrix 
elements H(P) have values =1 and -1), the element-by- 


element multiplication of the matrix Tr(X) by given 


matrix H(P) having the same dimension. These element- 
by-element multiplications are actually a change in the 


sign of the elements of the matrix Tr (X), for which the 


corresponding elements of the matrix H(P) have a value of 
-1. Thus, for the volume CTrs of computational operations 
when using Transpositions_matrix, we can record the 
following estimate 


CTrs = | norm + n? bitwise multiplications + n? changes of 
sign (11) 


Since bitwise multiplications and changes of sign of 
integer values takes less microprocessors cycles than 
cycles of floating point multiplication, the comparison of 
the volumes of the computational operations CH and CTrs 
shows that 


CTrs< CH (12) 
and therefore the method for creating a 4- and 8- 
dimensional orthogonal matrix R that rotates a given 
vector X in the direction of the x; coordinate axis that uses 


Transpositions_matrix has higher performance than the 
method that uses Householder Reflection 
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3.3 THE COMPUTATIONAL EFFICIENCY OF THE 
METHOD THAT USES BLOCK DIAGONAL 
TRANSPOSITIONS MATRICES 


In [1] the possibility of using Transpositions matrices to 
obtain an orthogonal matrix that rotates a given N- 
dimensional vector X in the direction of a given coordinate 
axis (i.e. axis xı) is considered. Block diagonal 
Transpositions matrices are used for this purpose. This 
paper examines the effectiveness of a method that uses a 
block-diagonal matrix Trb, consisting of two blocks of Trs 
matrices. When multiplying the block-diagonal matrix Trb 
by the specified vector X, it is obtained a vector Y1, which 
contains two non-zero elements. By multiplying the vector 
Y1 with a Givens matrix, is obtained the vector Y2, which 
is in the direction of a coordinate axis xı. Obtaining this 
result can be represented as follows: 


Y2=G(1,m +1,0).Y1 =G(1,m + 1,6).Trb.X 
Trs(X1/|X1) 0 
0 Trs(X2/|X2I|) | 


=G(1,m +1,0). 


=R.X =|]x\,0....0]" 


where: 
e X=[X1,X2]T is given vector of dimension N < 16. 


e X1, X2 are parts ot input vector X=[X1,X2] of 
dimensions m, n, N=m+n, m=2p, n=2q, p,q [2, 4, 
8] 

e = Trs(X1/|X1||) and Trs(X2/||X2]|) are Transpositions 
matrices [1] 

e G(1,m+1,6) is Givens matrix, which rotates vector 


Y1 in the direction of axis x1. 


3.3.1 ALGORITHM FOR CALCULATING V 
MATRIX USING BLOCK DIAGONAL Trg 
MATRIX: 


The algorithm for calculating V matrix using block 
diagonal Trb matrix, includes the following operations: 
Step 1: Split the given vector X into vectors X1 and X2 
with dimensions 2, 4 or 8. 

Step 2: From the given vector X = [X1, X2]" a vector x 
is obtained. The two parts of vector X - the vectors 
Xland X2 , are calculates by vectors X1 and X2 as it is 


given in formulas: 


X1 = X1/|X]| 
X2 = X2/|x2| 14) 
Step 3: The matrix V is obtained as follows: 
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V =G, (1,m+1,0).Trb 
Trs(X1) 0 


=G,(1,m+1,9). — 
mee | 0 mgo) (15) 


where: 


e Gn(1,m+1, 0) is Givens matrix, which has 
coefficients 


c=|xi]/|x 


S-KA ao 


e Trb is block diagonal matrix, which consists blocks 


Trs(X1) and Trs(X2) as it is given above. 


Multiplying the matrix V obtained by formula (15) by the 
vector X gives as a result vector Y=[|{X||,0,...,0]' . The 
matrix V calculated by formula (15) has the form: 


m x 


CX, CX, CIm | (“SX mr SIm 3X, 
Xa -X -Xa 0 0 0 
oe Ze Xei =X] 0 0 0 
~ | să% SX | 3 oe ee Gx. 
0 0 0 | Xm Ima e Xa 
L 0 0 0 Xa X pas Ean | 


Fig. 1 Scheme for computing the V matrix 


As can be seen from Fig. 1, obtaining the V matrix 
involves obtaining two Trs matrices of sizes m and n-m, 
respectively, and a total of 2n multiplications with the 
coefficients of C and S of the Givens matrix. In addition, 
to calculate the coefficients of C and S of the Givens 
matrix, it is necessary to calculate the norm ||X|| of the 
given vector X. Thus, for the volume Cr of computational 
operations when using block diagonal Transpositions 
matrices using formula (11), we can record the following 
estimate 


Cr = 2.(1 norm + n?/4 bitwise multiplications + n?/4 
changes of sign) +1 norm 


and consequently 


Cr = 3 norms+ n?/2 bitwise multiplications + n?/2 changes 
of sign (17) 


A comparison of estimates of the volume of calculations 
using Householder Reflection and using block diagonal 
Transpositions matrices shows that the volume of 
calculations is approximately the same Cr ~ Cy. In some 
cases, when using a more efficient vector norm calculation 
algorithm, using block diagonal Transpositions matrices 
may be a better choice. 


Theorem: V matrix, created using algorithm, described 
above, is matrix of rotation. 
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Proof: The V matrix is represented as a product of a 
Givens matrix and a Trb matrix V=G,(1,m+1,0).Trb (15). 
Since the Givens matrix is a rotation matrix, for the V 
matrix to be a matrix of rotation it suffices to prove that 
the Trb matrix is a matrix of rotation. For this purpose, it is 
necessary to prove that Trb is an orthogonal matrix (i.e. 
Trb.Trb™=I, where I is identity matrix) and that det(Trb)=1. 


First, for Trb.Trb™= we get: 


Trb.Trb! Ea 1 Ba 0 | 


Trs X2) 0 Trs(X2)" 
Z Ee 0 | 
0 Trs(X2)Trs(X2)" 


whence, taking of consideration that Transpositions 
matrices are orthogonal matrices we get Trb.Trb™=I. 


Second, since the matrix Trb is a block diagonal matrix 
whose blocks are Transpositions matrices Trs(X1) and 
Trs(X2), then for the determinant of the matrix Trb we 
get det(Trb) = det(Trs( X1)).det(Tr s(X2)), whence, 
taking into account that Transpositions matrices are 
matrices of reflection [1], ie. det(Trs(X1))=-1 and 
det(Trs( X2)) = -1 , we get det(Trb)=1, Q.E.D. 

The following example shows the code of a Matlab 
function that creates and returns orthogonal matrix which 


rotates a given vector X in the direction of the coordinate 
axis xı using block diagonal Transpositions (Trs) matrices. 


function V = OrthogonalTrs(X1,X2) 

nl = length(X1); n2 = length(X2); 

ifn2~=1 && nl ~= 2 && nl ~=4 &K& nl ~=8 
error(‘Dimension of X1 have to be 1,2,4 or 8’); 

end; 

if n2 ~= 1 && n2 ~=2 && n2 ~= 4 && n2 ~=8 
error('Dimension of X2 have to be 1,2,4 or 8°; 

end; 

end; 

X=[X1;X2]; %Full input vector 

normX=norm(X); %Norms of vectors 

norm1=norm(X1); norm2=norm(X2); 

Xin = X1/norm]; 

X2n = X2/ norm2; 

Trs1 = TrsMatrixs(X1n); %Transpositions matrices 

Trs2 = TrsMatrixs(X2n); 
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V = [Trs1 zeros(n1,n2); zeros(n2,n1) Trs2]; % Orthogonal 
Trs matrix 


C=norm1/normX; S=norm2/normX; %Coeficients of 
Givens rotation 


V(1,1:n1) = C* X1n'; 
V(nl+1,1:n1) = -S* X1n’; 

VC, nl+1:nl+n2) = S* X2n'; 
V(nl+1, nl+1:nl4+n2) = C* X2n'; 


Example 3 Matlab function, that creates orthogonal 
matrix, which rotates a given vector X in the direction of 
the coordinate axis x; (V matrix) using block diagonal 
Transpositions (Trs) matrices. 


As it can be seen, function OrthogonalTrs(X1,X2) 
calculates norms of input vectors X1 and X2, obtain V 
matrix as block diagonal matrix, the blocks of which are 
Transpositions (Trs) matrices. The rows number / and 
nl+1 of block diagonal matrix are replaced by the result of 
multiplication of vectors X1 and X2 by the coefficients C 
and S of Givens matrix, calculated as it is given in formula 
(16). 


IV. NUMERICAL EXPERIMENTS 


Executing the Matlab code below yields a V matrix using 
the OrthogonalTrs function. Obtained V matrix rotates 
given vector X=(1,2,3,4,5,6)' in the direction of coordinate 
axis xı. The vector X is split into vectors X1=(1,2,3,4)7 
and X=(5,6)" of dimensions 4 and 2. When executing the 
following Matlab code 


X=[12345 6]; 

X1=[1 2 3 4]'; X2=[5 6]'; 

V = OrthogonalTrs(X1,X2); 
Y=V*X; 


is obtained the following orthogonal V matrix: 


| A Array Editor V - o x 
| Fike Eda View Web Window Help 
Numeric format [shonG v| Size: fe by |6 x| 
Ja pr Cae Make a 


1f 0.10483) 0.20966| 0.31449 0.41931) 0.52414 0.62897 


2| 0.36515] -0.18257| 0.7303, -0.54772 0 0 
3 0.54772 ~0.7303| -0.18257 0.36515) o 0 
4 0.7303 0.54772| -0.36515 -0.18257 0 0 
s| -0,14948| -0,29896 -0.44844 -0.59792| 0.36757] 0.44109 


6 0 0 0 0| 0.76822) -0.64018 


Fig. 2 An orthogonal V matrix obtained by the 
OrthogonalTrs function. V matrix rotates the vector 
X=(1,2,3,4,5,6)" along the x; coordinate axis. 
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As a result of multiplying the obtained matrix V by the 
specified vector X, a vector Y is obtained, which is in the 
direction of the coordinate axis xı. 


V. CONCLUSION 


The proposed algorithm for creating an N-dimensional 
orthogonal matrix V, which rotates a given vector X of 
dimension N = m + n, m=2?, n=24, p,q € [2, 4, 8] in the 
direction of the x; coordinate axis using Transpositions 
matrices (Trs) allows to obtain similar performance 
compared to the computational efficiency of the method 
that uses Householder Reflection under the above vector 
dimension limits. The matrix, obtained by the proposed 
algorithm, can be useful in performing some operations in 
linear algebra such as QR decomposition of matrices. 
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